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We study the creation of vortex states in a trapped Bose-Einstein condensate by 
a rotating force. For a harmonic trapping potential the rotating force induces only 
a circular motion of the whole condensate around the trap center which does not 
depend on the interatomic interaction. For the creation of a pure vortex state it is 
necessary to confine the atoms in an anharmonic trapping potential. The efficiency of 
the creation can be greatly enhanced by a sinusodial variation of the force's angular 
velocity. We present analytical and numerical calculations for the case of a quartic 
trapping potential. The physical mechanism behind the requirement of an anharmonic 
trapping potential for the creation of pure vortex states is explained. 



I. INTRODUCTION 

The experimental realization of magnetically trapped Bose-Einstein condensates (BEC) of alkali- 
metal atoms has led to a broad interest in the properties of this state of quantum matter. One 
aspect of particular interest is to study the superfluid behaviour of these systems. An ideal tool to 
examine this is to rotate the condensate and to observe whether such a rotation can result in the 
creation of vortex states. 

Theoretically many aspects of vortex states in a trapped BEC have been discussed in the literature. 
The critical angular velocity required to create a vortex state in trapped BECs was derived by Baym 
and Pethick [Q. Dalfovo and Stringari determined numerically the shape of the vortex state ||^ 
and the moment of inertia in harmonic traps jj]. Ho and Shenoy j5| discussed the influence of the 
internal atomic structure on vortex states and Barenghi |^ studied vortex waves. Wilkin et al. ^ 
considered the fragmentation of a rotating BEC with attractive interaction. Dodd et al. Q have 
numerically studied excitations around a vortex state. Rokshar showed that the stability of a vortex 
state depends on the shape of the trap ||] and considered the properties of vortex states on a torus 

0- . . . 

Despite recent interest in the rotational properties of BECs, the experimental realization of a vortex 

state remains to be done In contrast to superfluid Helium, the trapped BECs are not in direct 
contact to an external container. How to rotate the trapped BECs into vortex states is still an 
open question. We have recently proposed to transfer angular momentum to a BEC by using multiple 
travelling wave laser beam s |12[ and to create vortex states by means of Laguerre-Gaussian laser beams 
using a Raman transition |]13[^^A related idea based on Raman transitions has also independently been 
proposed by Bolda and Walls iQ. In addition, a combination of multiple laser beams with Raman 
transition has recently been employed to discuss the creation of dark solitons and vortices in trapped 
BECs ||. 

In the present paper we study different aspects of the quantized circular motion of a trapped BEC 
driven by a spatially homogeneous rotating force. In Ref. we showed that such a force can be 
created by four travelling wave laser beams. In particular, we analyse the reason why in a harmonic 
trap the rotating force only leads to a rotational state of BECs which is a coherent superposition of 
vortex states. On this basis, we further demonstrate how to create a pure vortex state by introducing 
an anharmonic trapping potential. 

The paper is organised as follows. In Sec. II we derive the Gross-Pitaevskii equation for trapped 
BECs composed of two-level atoms in the presence of highly detuned laser flelds with different fre- 
quencies. In Sec. HI, we review the proposal of Ref. |jl^ and apply the general model to the transfer of 
angular momentum to BECs by a configuration of four travelling-wave laser beams. The force induced 
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by these laser beams is spatially homogeneous and rotates with the frequency difference between the 
laser beams. In Sec. IV we discuss the center-of-mass motion and the shape of a harmonically trapped 
BEC driven by the rotating force. Further we analyze the physical mechanism of why rotating a har- 
monically trapped BEC cannot create a pure vortex state. The possibility to create a pure vortex 
state of an ideal BEC is demonstrated in Sec. V by introducing an anharmonic trapping potential. 
In Sec. VI we consider an interacting Bose gas in an anharmonic trap. With a sinusodial modulation 
of the frequency difference between the laser beams the nonlinear energy shift due to the interatomic 
interaction can be compensated. We numerically show that an almost pure vortex state can be cre- 
ated. The conlusion is included in Sec. VII with a short discussion of the relation between the present 
proposal and that of Refs. [|l3|,E4j based on Raman transitions. 



II. GROSS-PITAEVSKII EQUATION FOR A BEC IN MULTIPLE LASER BEAMS 

We consider a condensate of N two-level atoms at zero temperature confined in a trapping potential. 
In the formalism of quantum field theory this system can be described by the Hamiltonian 

The symbol /i denotes the chemical potential. The center-of-mass Hamiltonian 

i?c.,n. ^-^^+V^trap(x) (2) 

contains the kinetic energy and the trapping potential Vtrap(^)j which is assumed to be cylindrically 
symmetric around the z-axis. '^i{x) denote the field operator for atoms in the internal ground state 
(i = g) and excited state [i = e) with internal energy Ei. M dentoes the atomic mass. Without 
affecting the generality of our final results we have assumed that the trap potential is the same for 
both internal states. The second integral in Eq. (|l|) describes the interatomic interaction in the 
absence of laser fields. The factors gij are defined as Airh^aij/M, with being the s-wave scattering 
lengths for atoms with different internal states. 

The interaction of the atoms with the laser beams is incorporated by the electric dipole coupling in 
the rotating wave approximation, 

/n 
d3x{vl/t(a-)vl/g(a-)^nr!(+)(a-)e-''^"* +i?.c.} . (3) 
a=l 

Here cja, o. — 1, . . . , n is the frequency of the ath laser beam, and ^^a'\x) is defined as d ■ E^a'\x)/h 
with d being the dipole moment of the atoms and -Ea^"* denoting the positive frequency part of the 
ath laser's electric field. 

It has been discussed at length in the literature |pj|-p^ how this Hamiltonian can be reduced to an 
effective Hamiltonian for atoms being in their internal ground state in the presence of a single laser 
beam. For this reason we will only briefly discuss how this can be done in the presence of several laser 
beams with different frequencies. In order to avoid spontaneous emission all laser beams have to be 
detuned far off resonance so that the population of the excited-state atoms is negligible. In this case 
the Heisenberg equation of motion for the excited-state field operator can be adiabatically eliminated 
by assuming that the detunings := LOa — {Ee — Eg)/h of the laser beams from the atomic resonance 
frequency are all in the same order of magnitude and are the largest characteristic frequencies of the 
system [|6j. The resulting expression for is given by 

(4) 

The Kernel L{x — y) describing the photon-exchange induced dipole-dipole interaction is defined in 
Ref. pq . The parameter 7 denotes the sponataneous emission rate from excited atoms and is assumed 
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to be much smaller than the detuning of the laser beams in this paper. The expression (Q) for can 
be inserted into the Heisenberg equation for '^g thus leading to the effective equation 

in4>g{x,t) = ~[Ho,'ifg{x,t)] + VL{x,t)'fg{x,t) + VNL{x,t)-fg{x,t) , (5) 
with the light-induced linear potential being given by 

h " 

^^^"^'^^ ■= 4(A + z7/2) ^ i^i^Hmi'\£)e-"^^-'^'^' ■ (6) 

Note that in the denominator we have replaced Aa by the average detuning A of all lasers. This has 
to be done for consistency since otherwise Vl would be non-Hermitian even in absence of spontaneous 
emission. Mathematically this can be made because the difference A;, — between two detunings 
has to be treated as of higher order in 1/A. Assuming that the condensate has a coherence length 
larger than the characteristic distance of the photon-exchange induced dipole-dipole interation, we 
derive the nonlinear potential 

VMS, t) = g^g^lix, t)^gix, t) 4(A^;,V4) 

The strength of this potential depends on the intensity of the applied laser fields. As we will be 
concerned only with very weak laser beams Vnl can safely be neglected in the present case. 

To describe the collective evolution of a condensed gas of bosonic atoms we will adopt the mean 
field approach in Eq. (j^). This approximation describes a trapped Bose-Einstein condensate quite 
well and amounts in replacing the field operator ^g{x,t) in Eq. (||), after the calculation of the 
commutator, by a collective wavefunction y/Nij}{x,t) of all atoms in the condensate, with ^/i being 
normalized to one. This procedure results in the Gross-Pitaevskii equation for ?/^, 

'^^^^^ = {H,.„,,+VL{S,t)+gN\4,{x,t)\^}4>{x,t) . (8) 

For brevity we have introduced g := ggg. This equation determines the nonlinear dynamics of trapped 
BECs in multiple laser beams and defines the effective Hamiltonian 

H ^H,,^,+VL+Ng\,l^\''i, (9) 

for the collective wavefunction tp. 



III. TRANSFER OF ANGULAR MOMENTUM BY AN OPTICAL POTENTIAL 



The design of an optical potential Vl which can rotate the BEG by transferring angular momentum 
was addressed in Ref. to which we refer for all details of the construction. Here we will only briefly 
sketch the important features of this optical potential. 

The potential Vl {x, t) of Ref. jl^ is created by four traveling wave laser beams of equal intensity 
and width much larger than the size of the condensate (this is usually the case) so that they can be 
considered as plane waves. The laser beams 1 and 2 propagate along the x- and y-direction have the 
same frequency loi = Ld2 = lo- Their phases are chosen to be equal at the origin so that the Rabi 
frequencies of the two lasers have the form Vl']^\x) = Vt^eyi^likx] and Vl^^\x) = rioexp[zfcy]. The 
other two laser beams have the frequency W3 = (jJ4 = lo' and propagate in directions slightly different 
from the x- and y- direction with a small angle 9. The frequency difference w — w' is assumed to be 
so small that the difference between \ka\ and \kb\ for any a, 6, = 1, . . . , 4 is negligible. The wavevector 
of the third laser beam is given by kcos{9)ex — fcsin(6')ey ~ ke^ — Skcy, where 

Sk := kO (10) 

denotes the deviation of the wavevectors. Similiarly the 4th laser's wavevector is given by key — 
Sktx- The phases are chosen so that laser 4 is in phase with lasers 1 and 2 at the origin whereas 
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the phase of laser 3 is shifted by — 7r/2. The third Rabi frequency is then given by fig ' {x) = 

exp[— j7r/2]f2o exp[i(fca; — Sky)] and the fourth by fl^^"^ (x) — f2o exp[j(fcy — Skx)]. 
The described configuration of laser beams leads to the optical potential 

where we have introduced the effective Rabi frequency flcff '■— |rioP/(4A) and the time dependent 
phase (j)(t) := t{uj — uj'). In this section and in the rest of the paper we neglect c-number terms in the 
optical potential since they only produce an overall phase shift in the wavefunction. 

Since the wavelength of an optical laser is much smaller than the size of a typical condensate (several 
micrometers) all exponentials in Eq. ( pT| ) containing the wavevector k are rapidly varying over the 
condensate and therefore average to zero The remaining exponentials, exp[ii5fca;] and exp[i5A;?/], 
can be approximated by expressions linear in 5k if the angle 9 is small enough. The resulting optical 
potential then takes the form. 

Vl{,x, t) = 2nrJeff'5fc{sin((/)(t))x - cos((j){t))y} . (12) 

It obviously produces a rotating spatially homogeneous force over the size of the condensate. The 
rotating force can transfer angular momentum to the condensate. This can be seen easily by rewriting 
it in the form 

Vl{x, t) = inVl^s5k{{x + iy)e-*"*(*) - {x - iy)e"l''^^^ . (13) 

Because of the commutation relation [L^, {x ± iy)] = ±h{x ± iy) the operator x -\- iy increases the 
orbital angular momentum by Ti. Thus, if the frequency difference w — w' is nearly resonant with a 
transition between two trap eigenstates, the angular momentum of the condensate is increased by the 
force @. 

Since Vl is produced by laser beams we have to estimate the rate of decoherence of the condensate 
due to the spontaneous emission of photons which is given by the probability that an atom is excited 
times the spontaneous emission rate 7. The excitation probability for off-resonant light is of the order 
of |r2op/A^ = Oeff/A. The decoherence time is therefore given by A/(7r2cff) and is therefore much 
larger than 1/ricff since the detuning is much larger than the natural linewidth. For instance, if \/^lcs) 
is about 10 ms then for A = IOO7 the decoherence time is in the order of a second. 

One also has to address the question of how stable the frequencies of the laser beams must be so that 
ui — ui' is well defined. It turns out that this poses no problem at all since only the frequency difference 
Lo — uj' \s required to be stable. Since in the proposal of Ref. |l2j all four laser beams originate from 
the same laser source, and since the frequency difference is created by modulating the laser beams, 
it is clear that only the modulator frequency must be stable. The stability requirement on the laser 
beams can therefore easily be fulfilled by any common laser source. 

We also point out that laser beams are of course only one possibility to produce a rotating homo- 
geneous force. The same potential could be created by spatially varying and rotating magnetic fields, 
for instance. The following analysis of the condensate's motion in the potential (|l^) therefore covers 
this case as well. 



IV. HARMONICALLY TRAPPED CONDENSATE: COHERENT ROTATION 

We now turn to the effect of the optical potential (|l^ on a condensate trapped in a harmonical 
potential of the form 

V,..,{x-)^^^lz' + ^ul{x^ + y'). (14) 

ujz and wj, denote the trap frequency in the z-direction and in the x-y-plane, respectively. 

We have shown in Ref. that in the case of an ideal Bose gas and a resonant optical potential 
(i.e., UJ — uj' — uj\_) the condensate rotates around the center of the trap with increasing mean radius 
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but preserving its shape (comp. Fig. 1). It was also demonstrated that this motion corresponds 
to a coherent superposition of vortex states ipn (oc exp[m(/3], where ip is the azimuthal angle in the 
x-y-plane) with different winding number n. 

In this section we extend this result to an interacting Bose gas. We start from Eq. (||) including 
the resonant optical potential ([1^). As is well-known the time evolution of the expectation value of 
any operator O is governed by the equation 

^(0)4([i^,0]). (15) 

It is straightforward to derive from this relation the corresponding equations for the position operators. 
Observing that / d'^xji/'p^ili/'P, i = x,y,z vanishes the center-of-mass motion takes the simple form 

f^(0+u;l{0^m, (16) 

for = x,y,z), where the hght-induced force f{t) is given by f{t) = 2{hilcS^k/M){— sm{Lj±t)ex + 
cos{uj±t)ey}. For a BEG in the ground state the condensate's mean position and mean momentum 
are initially zero ((.t)(0) = , {p){0) ~ 0) so that wc find for the solution 

(x) = RQ{ex[uj±_tcos{uj±t) — sin(wj_i)] + CyLu ±t sin^uj ±t)} , (17) 

withi?o := hn^gSk/iMul). A gain the condensate performs a rotation around the trap center whereby 
the radius of the rotation increases linearly in time. This result was to be expected since usual 
interatomic potentials are translation-invariant and therefore affect the center-of-mass motion only 
indirectly. It provides an extension of Ehrenfest's theorem for interacting Bose condensates. We 
remark that the result ( pT| ) deviates from the corresponding result of Ref. by the factor — sin(ti;_Li) 
inside the square brackets. This difference arises because in Ref. fl^ a rotating wave approximation 
(RWA) has been made. As usual for RWA the difference is negligible for large times uj±t ^ 1. 

To demonstrate that, as in the noninteracting case, the shape of the condensate is not affected by 
the optical potential we consider the time evolution of the variance of the position operators. Using 
J d^x^\i}\'^d^\ip\'^ = - Jd^x\^\^/2 and defining A^ (^p^ + PiO - 2(0 (P?) we arrive at 

dt M 

^ = ^{Ap, f - 2Mul{Mf + g J d'umu)\* (18) 
d{Ap^)^ 



dt 



MujiA^-ihg / d^u{(a^V*(w)) (w) - ^^-c.} 



Because neither the center-of-mass motion nor the optical potential enter into these equations we can 
infer that the shape of the condensate remains unaffected by the optical potential. 

It is worth noting that one can derive from Eq. ( p^ ) the mean energy conservation law 

^(Ap)2 + ^MiculiiAxf + (Aj/)2] +c.2(Az)2} + g J d\\ij{u)\^ = const. , (19) 

which nicely demonstrates the relation between the interaction energy and the widths of the conden- 
sate's wavefunction. 

We conclude that, albeit the potential ([l3| ) seems to be perfectly suited for this task, neither a pure 
vortex state nor even a vortex line far away from the center of the trap is created. However, this does 
not rule out the possibility to form a vortex state during the formation of the condensate by directly 
cooling the gas in a rotating trapping potential. Our result may also be helpful to understand why 
the recent experiments at MIT |ll|] attempting to observe persistent flows around a potential wall at 
the trap center by rotating the trap cannot achieve the expected result. 

To understand why the potential ( p^ ) produces a superposition of vortex states rather than a pure 
one we observe that it couples only neighboring eigenstates of the trap. This can be seen by expressing 
the coordinates x, y in terms of the corresponding annihilation operators of the trap (i.e., x (x Ux + al., 
for instance). Since the potential is linear in these operators it couples states with energy nhLu± (in 
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the noninteracting case) to states with energy {n + l)hLLj±. If the coupHng is resonant {4>{t) = uj±t) 
the potential not only resonantly transforms the ground state to the first excited vortex state, 
but it also resonantly transfers the nth vortex state (with energy nhuj±) to the {n+ l)th vortex state, 
because of the equal energy interval between these energy eigenstates. As a result, the long-time 
evolution of the system driven by the rotating force will lead to a superposition of vortex states with 
different winding number. 

This implies that to create a pure vortex state one must destroy the equal interval between the trap 
energy levels to avoid the simultaneous excitations of different vortex states. We will analyze this in 
the next section 



V. CREATION OF A VORTEX STATE IN AN IDEAL BOSE GAS 

In the preceding sections we have presented a scheme to transfer a harmonically trapped condensate 
into a rotating state which is a coherent superposition of vortex states. We also analysed why the 
scheme does not create a pure vortex state if the trapping potential is harmonic: the equal interval 
between the energy levels in a harmonic trap causes the condensate to climb up a ladder of resonant 
transitions into different vortex states. 

The situation is different if an anharmonic (yet still rotationally symmetric) potential I4rap is used 
in Eq. (||) for which the interval between the energy levels is not equal. Therefore, if the frequency 
a; — w' in the optical potential (|l^ ) is choosen to be in resonance with the transition from the ground 
state to the first excited vortex state, the transition to the next vortex state will in general be out of 
resonance (see Fig. 2). This fact can be exploited to prepare a pure vortex state in an anharmonic 
trap. 

To estimate this quantitatively we first consider the ideal Bose gas and expand the wavevector in 
terms of the simultaneous eigenstates ipm,s of the Hamiltonian (|^) and the orbital angular momentum 
Jz, i> = I]m,sexp[-zTO(/)(t)]c„,s(i)'/'m,s- The eigenstates fulfiU the relations Hc.m.i>m,s = E^,si>m,s 
and Jzipm,s = 'mhtpm,s- It should be noted that the inclusion of the phase factor exp[— zm</>(t)] in this 
expansion is equivalent to the description of the system in a frame of reference rotating with angular 
velocity 4>(t) = to — lu' (where we left open the possibility that the frequencies can vary in time). 
In this description the broken symmetry between states of positive (m > 0) and negative (m < 0) 
angular momentum becomes obvious. Since the latter are counterrotating for uj — uj' > transitions 
to these states are highly suppressed. 

Using this expansion the Schrodinger equation (H) for the ideal condensate with g = can be 
reduced to 

ihCjn,s = {Ejn,s - mh4l(t)}Cm.s + fl^cS^k ^ {gm;s,s' Cm- + 9m+l:s',sCm+l,s' } (20) 

s' 

With qm-s,s' ■= J (Pxil)*^ ,.{x + iy)ipm-i,s' ■ 

We assume that, as it is the case for the harmonic potential. Hem. has a unique ground state with 
energy i?o,o and that the energy levels i?i,s are well separated, with Ei^q being the lowest vortex 
excitation. In this case, if we choose oj — oj' = (-Bi,o — Eo^o)/h the transition from V'0,0 to -01,0 is in 
resonance while the transition from ipi^Q to 1P2.S is detuned by the amount As := uj — uj'—{E2.s — Ei^o) = 
-£■2,8 + £-0,0 ~ 2i?i_o- Furthermore, if the condition |f^cff<^fc'72;s,o| is fulfilled, the transition from 

■01,0 to 7/12.5 will be highly supressed so that we effectively have a two-level system consisting of the 
states "00,0 and ipifl. It is obvious that after a time -k/v the complete population will have been 
transferred to the first excited vortex state, where we have defined the Rabi frequency 

v := |ilcff'5fcgi;o,o| (21) 

of the transition between -00.0 and 4'i,o This situation is mathematically very similar to the resonant 
excitation of the first excited state in an atom by a laser beam. In this case one can safely neglect all 
states beside the ground and the first excited state, too. 

As a specific example we consider a two-dimensional anharmonic oscillator for which the trapping 
potential is given by Vtrap — i^f'^ , where r is the two-dimensional radial variable and k determines the 
strength of the potential. The neglection of the z-direction does not alter the predictions for the ideal 
Bose gas and is justified for the interacting Bose gas if excitations along the z-axis are sufficiently 
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supressed. Approximate eigenstates of the resulting Hamiltonian can be obtained by application of 
Ritz's variational method. For later use we include the nonlinearity and minimize the expression 



E 



(22) 



for the conserved energy. The coupling constant giD that we use for an effective description of 
the vortex creation in two dimensions is related to the three-dimensional coupling parameter g by 
92D = gN/li^. Using the ansatz 



n+1 



7rn! 



72 



(23) 



we find for the variational energy eigenvalues the expression E, 



'n,0 



E. 



(NL) 



with 



E. 



n,0 



.(n + l) 



(3 + Gfn 



2{n + 2) 

1 + Gfn 



1/3 



E. 



(NL) 
n,0 



£Gfn{n + l) 



2(n + 2) 

1 + Gfn 



1/3 



(24) 



e!^q represents the contribution of the linear part of the Hamiltonian (kinetic and potential energy) 

and E^f^'^ denoted the contribution of the nonlinear interaction energy. We have introduced the 

energy scale £ := kZq, the dimensionless interaction parameter G '■= -^ff2D/(27r?i^), and the numerical 
factor /„ :— {2n — l)!!/[(n + 1)!2"]. The characteristic length scale of the system is given by lo := 
[a^/(2M/c)]i/^ For the ideal Bose gas we set G ~ and find for the lowest energy eigenvalues the 



values E: 



0,0 



2.38 Eio^ 5.45 E; 



2,0 



a 9- £, and Eq i sa 9.64 • £. We have determined the energy 
Jqj. ^j^q wavefunction. The largeness of Eq i allows us 



£^0,1 by using the ansatz -0o 1 oc (e — r )e' 
to neglect all transitions to this state. 

The interaction matrix elements are of the order of Iq for these four states, but we will need only 
the value for (?i;o,o ~ 0.79 • Iq and q2;o,o ~ 1-05 • lo- The transition from |0, 0) to |1,0) will be in 
resonance if we choose the frequency difference <j) = 3.07 • £/h. The transition from |1,0) to |2,0) 
is then supressed if the interaction matrix element hflcffSk ■ g2;0,o is much smaller than the energy 
difference £'2,0 ~ -^'i.o — h4> ~ 0.48 • £. This is guaranteed if we set the transition energy -Etrans '■= 
equal to £ /2Q. The laser scheme presented above then will create a vortex state after a time tt/v. We 
have depicted this physical situation (to scale) in Fig. 3. 

To estimate its experimental realizability and to compare it with the interacting Bose gas we consider 
the case that the numerical value of the trapping potential strength n is given by 7.66 • 10~* J/m'* 
so that for ^''Rb atoms (M = 1.45 • 10~^^ kg) the trap size turns out to be /q = 8.91 • 10~'^m (These 
numerical values turn out to be convenient for the numerical calculations presented below). We then 
find for the energy scale £ « 4.83 • 10~^^ J and for the characteristic time for vortex creation -k/v k, 70 
ms. This well within the range of current experiments with Bose-Einstein condensates. We mention 
that we have also made a numerical study of this creation process for the numbers given here. The 
result is in excellent agreement with our analytical calculations but will not be presented here since 
the analytical results already perfectly describe the physical situation for a non-interacting Bose gas. 
For all details of the numerical procedure we refer to the next section where it is applied to the 
correponding case of an interacting Bose gas. 



VI. CREATION OF A VORTEX STATE IN AN INTERACTING BOSE GAS 

To study the infiuence of the atomic interaction analytically we assume that the nonlinearity in Eq. 
(jsj) is weak enough so that we still can use the variational ansatz (^). This is not allowed for the 
strong interaction regime where this ansatz becomes meaningless, but it includes the case that the 
interaction energy is of the order of the kinetic and potential energy and thus overcomes the restriction 
of perturbation theory. 

We further assume that radial excitations of the condensate are unimportant because the en- 
ergy levels E„i^s {s > 0) are considerably larger than the corresponding lowest value Em,o for the 
same angular momentum mh. We then only need to take into account the states V'm.o- The 
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numerical results presented below further justify this assumption. Writing the wavefunction as 
ip = '^^Cn(t)ex-p[—inAujt]ipn,o{^) and projecting Eq. (||) onto ipm.o we arrive at 

ihCm = {-E^Q - mh,(j)}Cm + hi}csSk{qm-0,OCrn-l + 9m+l;0,oCm+l } 

+2£;i~o^)|c„pc„ + ^ <i|cfepc™ + E E 4QC„,+,_,<g . (25) 

Here we have introduced the energies eI^\ = i?^"^ := '2g2D J d'^x\^pm,o\'^\i^k,o\'^ for scattering between 

atoms in state k and to, and E^.^\ :— g2D J dl^xip!^ ^Tpl Qipifiipm+k-ifi for anomalous scattering of 
atoms. 

To solve Eq. ( p5| ) analytically we first make the assumption that only the two states V'0,0 ^-nd ■01,0 
are populated. Though this is not justified for a quantitative analysis it provides very useful physical 
insight into the effectiveness of the process of vortex creation. Eq. (|2^) then reduces to 

thco = {e^om + 2E^oT^\cof + 4a |ci|'} CO + hfl^sSkqlofiCi (26) 
thci = (sfg) + 2e['1''^\c,\^ + 4a |col' - H} ci + nn,sSkq,,o^oCo . (27) 

This form of the nonlinear Schrodinger equation is particularly suited to study the effect of a time- 
dependent variation of the laser frequencies 4> = lj — to' . It was originally proposed in Ref. to use 
such a variation to compensate for the nonlinear interaction energy between the atoms. Here we will 
go into more detail and examine the question of which variation is the most effective one to create a 
vortex state. 

Let us first assume that cj) is constant and chosen so that initially (|coP = 1 , |cip = 0) the 
transition from -00. to ipi^o is on resonance, i.e., 

= <o^-<o^+i^i:^o^-24r- (28) 

After some time the population of the two states will be altered (|coP < 1 , |cip > ) so that the 
magnitude of the terms in curly brackets in Eqs. ( p6| , ^T\ ) is changed. Thus, the 2-level system is 
driven out of resonance by the nonlinear atomic interaction and a complete transfer of the condensate 
from 00,0 to 0i^o becomes impossible. 

It is not hard to see that a time-dependent variation of (f) can compensate for this interaction- 
induced shift of the resonance frequency. In Ref. Q we proposed to vary w — w' linearly in time. 
While this works well if the rate of change is suitably chosen, it has the disadvantage that one needs a 
relatively long time to create a vortex state. This is similar to a Landau-Zener transition [^,0 which 
is destroyed by nonadiabatic transitions if the frequency's rate of change is too high (sec Ref. |22[ , 
for instance). To improve this situation we assume that is implicitely chosen so that the transition 
between -00,0 and 0i,o is always on resonance, 

= 4^„' - + i?foH|coP - IciH + 24;^^)|cip - 24;^^)|cop . (29) 
Under this condition the solution of Eqs. (|2^, |2^) is given by 

co(i) = e-^'C) cos{iyt) 

ci(t) = -!^e-^Wsin(z.t), (30) 

l9l;0,o| 

where the phase factor 9{t) is easy to calculate but not needed here. This solution corresponds to 
a resonant Rabi Oscillation between the two states. Hence the two-level model predicts that the 
creation of a vortex state can be as effective as in the linear case provided the nonlinear interaction 
energy is compensated for. The corresponding variation of the frequency difference can be found 
by reinserting Eq. (|3^) into Eq. ( p9| ) and turns out to be a simple sinusodial variation with frequency 
2v. Such a sinusodial variation of the frequency difference can be realized by simple frequency 
modulation techniques. 

The above analysis was made under the oversimplified assumption that we only need two levels to 
take into account. The influence of other modes can be studied in several ways. For instance, they 
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can be treated as an elementary excitation around the solution ( p(\ j . For the sake of brevity we will 
not discuss the analytical treatment here and present instead numerical solutions of the full nonlinear 
Schrodinger equation (||). 

Numerical solution: Our numerical procedure started from the ground state of the interacting 
condensate in the quartic anharmonic potential. This state was found by propagating a Gaus- 
sian trial wavefunction in imaginary time and normalizing it between each step. The formation 
of a vortex state was studied by using the modified split-step method presented in Ref. |23| . 
This consists in approximating the time evolution operator exp[—iAt{T + V)/h] for small At as 
ex-p[—iAtV)/{2h)]exp[—iAtTh]exp[—iAtV/{2h)], where T and V denote the kinetic and potential 
part of the Hamiltonian, respectively. The terms containing V are diagonal in position space and 
can easily be applied to the wavefunction. The term involving T is computed by transforming the 
wavefunction to momentum space by using the fast-Fourier-transformation algorithm (FFT), applying 
the operator, and performing the inverse FFT. 

The physical parameters for the trap and the Rb atoms are the same as given at the end of the 
previous section (trap strength k — 7.66- 10~^ J/m'', trap size Iq — 8.91-10"^ m, mass M — 1.45- 10~^^ 
kg). In addition, the two-dimensional nonlinearity parameter g2D was assumed to be 2.6 • 10"^^ J m^. 
For Rb atoms the scattering length is about 5 nm. Therefore, if we relate an effective 3-dimensional 
peak density p^D to the peak of the 2-dimensional wavefunction 1^20 by \ip2D\'^g2D = Psog we can 
infer that this choice of 020 corresponds to a peak density p^n of about 2 • 10^'^ cm"'^. The strength 
of the optical potential ( |l3| ) was determined by setting hD,csSk = 2 • 10^^^ J/m. It turned out that for 
this system a grid size of 64 x 64 points was accurate enough for the numerical study of the creation 
process. 

The most laborious part of the numerical simulation was to find the proper variation of the frequency 
difference 4>- It is not difficult to find an approximation of the physical parameters in Eq. (^9|) by use 
of the variational wavefunctions (p3|). However, in order to get satisfactory results, these parameters 
must be determined more precise than the variationally calculated values which deviate from the exact 
values by up to 10%. We finally found that the best results are produced if the sinusodial variation is 
delayed for a time to — 50 ms to compensate for effects of other modes than tpo,o and '01, o- Thus the 
frequency difference was given by 

'^max t < to , , 

^{(^max + Winin) + (t^max - t^min) C0s(^(t - io))} t > to ' ^ ^ 

with Wmax = 1430 s~^, Wmin = 1205 s~^, and p = 19.5 s"-^. 

The numerical results for this choice of parameters are shown in Fig. 4 and 5. Figs. 4a-d show the 
modulus squared of the wavefunction at different times of the evolution (0 ms, 87 ms, 142 ms, 247 ms). 
It shows that for a moderately interacting Bose gas in an anharmonic trap an almost pure vortex state 
can be created. Fig. 5 displays the phase of the wavefunction at the same times (the coarse-grained 
structure at the border of the figures represents numerical noise). Fig. 5b demonstrates that during 
the evolution two vortex lines are created, one of which finally becomes the vortex state in Fig.4d and 
the other slips out of the trap. This situation can be described reasonably well by assuming that the 
wavefunction is a superposition of the form ^ = co4'o,o + ciipi^ + C2ip2,0: where V'2,0 describes a vortex 
state with angular momentum 2h. Fig. 4b shows the corresponding deformation of the wavefunction. 
The modulus of the expansion coefficients Ci can be approximately determined by looking at Fig. 5 
since the closer one (two) vortex lines are to the origin the larger the value of ci (02) is. We also 
mention that the vortex lines rotate around the origin during the creation process. This can be seen 
in an animation of the creation process |2^ ]. 



VII. CONCLUSION 



We have presented a scheme based on four laser beams which induces a spatially homogeneous force 
to transfer angular momentum to a trapped condensate. For a harmonic trap this scheme is equivalent 
to a rotation of the trapping potential around its original position. In this case the condensate will 
be prepared into a superposition of pure vortex states, corresponding to a circular motion of the 
condensate as a whole around the original trap center. The reason for the creation of a superposition 
rather than a pure vortex state is that the energy eigenvalues of the trap are equidistant for a harmonic 
trap. Thus, transitions to neighboring states are simultaneously resonant. In an anharmonic trap this 
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phenomenon does not occur. If the intensity of the laser beam is sufficiently low one can tune the 
transition from the ground state to the first excited vortex state into resonance while leaving the 
transition to the second excited vortex state far off-resonant. The interaction- induced energy shift 
appearing in a non-ideal BEC can be compensated by a sinusodial variation of the rotation frequency 
of the force. It then is possible to create a pure vortex state in an anharmonic trap. 

Beside employing an anharmonic trap one alternative method to create a vortex state is to use 
Raman transitions induced by Laguerre-Gaussian laser beams, as it was propsed in Refs. p3| , p^ . 
In this case a single Laguerre-Gaussian laser beam already carries an orbital angular momentum 
which can be transferred to the condensate. It is of interest to understand why the proposal with 
Raman transitions also works for a harmonic trapping potential and is not restricted to anharmonic 
potentials. In the notation used in the present paper the basic mechanism of the scheme based on 
Raman-transitions is that laser beams of the Laguerre-Gaussian type couple the condensate, which 
is initially in the ground state |0,0) of the trap and the internal state |— ), to a vortex state |1,0) 
and the internal state |-|-). Here the state vectors are defined by tpm.s{x) — {x\m, s). The interaction 
Hamiltonian is then roughly of the form Vl oc |1,0) (g) \+){-\ (g) (0,6| + |0,0) (g) |-)(+| (g (1,0|. This 
combination of internal transition and excitation of the center-of-mass motion guarantees that the 
condensate never can occupy higher vortex states like |2,0) g) |±) since the condensate can only 
oscillate between the states |0,0) g) |— ) and |1,0) g) |-|-) in the Raman coupling. Hence selective 
excitations of the internal state by Raman transitions avoid the resonant excitation of the ladder of 
transitions to higher vortex states. 
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FIG. 2. This diagram shows the ladder of transitions between different trap energy levels due to the optical 
potential: a) For a harmonic trap all transitions are resonant, b) In an anharmonic trap only the first transition 
is resonant so that only the first excited vortex state is populated. 
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Eo,o 

FIG. 3. The lowest lying eigenstates for the anharmonic trap. If the frequency difference uj — cj' is tuned into 
resonance with the lowest transition and the optical potential (^3|) (represented by Straus) is weak transitions 
to higher states are highly suppressed. 




FIG. 4. Numerical simulation of the condensate's time evolution during the vortex creation. The modulus 
squared of the collective wavefunction is shown at different times (a: ms, b: 87 ms, c: 142 ms, d: 247 ms). 
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FIG. 5. The same as Fig. 4. The phase of the wavefunction is shown, varying between tt (white) and — tt 
(black). A vortex Une corresponds to the end of a sharp borderline between white and black. 
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